Sound Velocity and Meissner Effect in Light-heavy Fermion Pairing Systems 
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In the frame of a four fermion interaction theory, we investigated the collective excitation in 
light-heavy fermion pairing systems. When the two species of fermions possess different masses and 
chemical potentials but keep the same Fermi surface, we found that the sound velocity in superfluids 
and the inverse penetration depth in superconductors have the same mass ratio dependence as the 
ratio of the transition temperature to the zero temperature gap. 
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The study on Cooper pairing between fermions with 
different masses promoted great interest both theoreti- 
cally and experimentally in recent years. A new pair- 
ing phenomenon, the breached pairing or interior gap, 
was proposed in the study of light-heavy fermion pair- 
ing systems 0,0,0]. The light-heavy fermion pairing can 
exist in an electron system where the electrons are from 
different bands, a cold fermionic atom gas where a mix- 
ture of 6 Li and 40 K can be realized, and a color super- 
conductor with strange quarks. 

In this Letter, we focus on the Cooper pairing between 
two species of fermions with equal Fermi surfaces but 
unequal masses. Such constraint on the fermion pairing 
can be realized, for instance, in the mixed atom gas of 6 Li 
and 40 K by adjusting the numbers of the two species to 
be the same. It is recently found that, for such a system 
the ratio between the transition temperature T c and the 
zero temperature gap Ao depends only on the mass ratio 
a = m,b/m a between the two masses m a and m&Q, 



A 7T \l + a) 



(1) 



with the Euler constant 7. Without loss of generality, we 
assume mj > m a in the following. 

What is the behavior of the collective excitation in 
such a system, and does the simple a dependence still 
hold when we go beyond the mean field approximation? 
In a superfluid composed of neutral fermions, the low en- 
ergy excitation is the Goldstone mode which is directly 
related to the specific heat at low temperature. In a 
superconductor composed of charged fermions, the col- 
lective mode is associated with the Meissner effect. 

We consider a system composed of two species of 
fermions with attractive interaction, described by the La- 
grangian density in Euclidean space x = (r = it, x) , 
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where ipi(x) describe the fermion fields, g is the coupling 
constant, and fj, a and jib are the chemical potentials of 
the two species. We will take the units c = h = fcs = 1 
throughout the paper. 

We can perform an exact Hubbard-Stratonovich trans- 
formation to introduce an auxiliary boson field <p(x) and 



its complex conjugate <p*(x). By defining fermion field 
vector <F* = ( rp*, ipb ) in the Nambu-Gorkov space, the 
partition function of the system can be written as 



Z = / [d**][d*][d<£][d0*]e 



with the kernel K defined as 
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where (3 is the inverse of the temperature, (3 = l/T. 

When the interaction is turned off, the chemical poten- 
tials \i a and fj,b are regarded as the corresponding Fermi 
energies. Since we focus on the Cooper pairing with equal 
Fermi surfaces, the Fermi momenta of the two species are 
the same, pp = p b F — pf = ^/2m^JI^ = y/2mb(J>bi an d the 
number densities are also the same, n a = rib = n/2 = 
P f /(6tt 2 ). In this case, the breached pairing state which 
is unstable to phase separation 0] and LOFF stated is 
ruled out. 

For the zero range interaction in J5J), we need a regu- 
larization scheme. For a solid, we can suppose that the 
attractive interaction is restricted in a narrow momentum 
region around the Fermi surface, pf — A < |p| < pp + A 
with A <C pp, like the classical BCS theory and the study 
in where A serves as a natural ultraviolet cutoff in the 
theory. For a dilute fermionic atom gas, we can replace 
the bare coupling g by the low energy limit of the two- 
body scattering matrix^], namely 



m 
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with the s-wave scattering length a s , e, p = p 2 /(2m), and 
the reduced mass m = 2m a mb/(m a + nib). While the 
results from the regularization scheme with the cutoff 
A and the scheme with the scattering length a s are the 
same in weak interaction limit, the latter is also a good 
approach to study the BCS-BEC crossover where the in- 
teraction is strong. Without loss of generality, we adopt 
the latter in this Letter. 

In mean field approximation, the boson field (j) is re- 
placed by its vacuum expectation value A which can be 
chosen to be real, and the gap equation which determines 
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FIG. 1: A schematic description of the dispersion relations for 
the two branches of free fermions (dashed lines) and quasipar- 
ticles (solid lines), pp is the common Fermi momentum. 



A is derived from the minimum of the thermodynamical 
potential Q = —TliiZ. At zero temperature, it reads 



Air a* 
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where E p is defined as E p = (e p — /i) 2 + A 2 with the 
average chemical potential /i = (p a + /Zb)/2. In the weak 
coupling limit, the solution of the gap equation at zero 
temperature can be analytically expressed as 
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The fermion propagator Q = K 1 in the Nambu-Gorkov 
space can be explicitly written as 

re \ ( iVn ~^ P ) -nA 
gK " PH ^ n -^?-El ' (8) 

where i/„ is the fermion Matsubara frequency, are the 
Pauli matrices, and (; p are defined as 



(9) 



with 

,,2 



the free fermion 



relations £^ 



dispersion 

p 2 /(2m a ) - p a and ^ = p 2 / (2m b ) - The disper- 
sion relations of the quasiparticles can be read from the 
poles of the propagator. A schematic description for the 
two branches of quasiparticles is illustrated in Fig^ Due 
to the constraint of equal Fermi surfaces on the pairing, 
all fermionic excitations are gapped, as in the standard 
BCS theory with equal fermion masses. 

The transition temperature T c is determined by the 
gap equation at finite temperature and at A = 0, 
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where f(x) = 1/ (e x ' T + l) is the Fermi-Dirac distribu- 
tion function. In the case of weak coupling, the num- 
ber equations can be approximated by the condition 
m a p a = m b fi b which gives = 2(e p - p)/(l + a) and 



£p = a£p- With the standard trickQ, we can reobtain 
the a-dependence of the ratio T c /A as shown in Eq. 

For a = 1, we recover the standard BCS result 
T c ~ 0.57Ao- If n and a s are fixed, we have the relation 
between the two critical temperatures, 



T c (m a ,m b ) = y/m a /m b T c (m a ,m a ), 



(11) 



which means that the critical temperature for the mixed 
6 Li and 40 K system is about \A/7 of the one for the pure 
6 Li system. 

We now start to investigate the low energy collective 
excitation in the superfluid state. The spontaneous sym- 
metry breaking and the associated Goldstone mode in 
the superfluid state provides an effective field theory ap- 
proach for the collective excitationjsl Ifll . In our case, 
the particle number is conserved which corresponds to a 
global U(l) symmetry of the phase transformation, 



(12) 



with an arbitrary and constant phase cp. The nonzero 
condensate A of Cooper pairs in the superfluid state 
spontaneously breaks the U(l) symmetry, and corre- 
spondingly, a Goldstone mode which possesses linear dis- 
persion at low energy is expected to appear. The low en- 
ergy dynamics of the system at low temperature is then 
dominated by the Goldstone mode. Since all fermions are 
gapped in our case, the fluctuation of the amplitude of 
the order parameter at low temperature can be neglected. 
Therefore, we can write the order parameter field as 

(f>(x) = Ae 2l9{ - X \ (13) 

Taking the standard approach [9I llOj. we can transform 
nonperturbatively the fermion fields as follows, 

= u*v e{x \ = i>*(x)e- ie w. (14) 

The transformation is designed to eliminate the phase 
fluctuation 9{x) dependence from the off-diagonal pairing 
potential terms (?/>*?/>^0 -I- c.c.) in the Lagrangian. The 
9{x) dependence can appear only in the kinematic terms 
of the fermion sector, and the Lagrangian density of the 
system can be expressed as 

. / D 2 \ ~ - - A 2 

C = xJj* -D T + lM\i>i+ {&4>b4>a + c.c.) + — 

y 2m, J g 

(15) 

with = dfj, + id p 9(x). Using the Nambu-Gorkov vec- 
tor, the partition function can be rewritten as 



Z 



'}[dnde]c J " dTld3 <^ Km+ ^ (16) 



with the kernel JC[9] defined as 

K\6] = K[A] - T 3 id T 9 - E_(V6») 2 + E+iV6» • V, (17) 
where the matrices £± are defined as 



l/m a 
±l/m b 



(18) 
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Since the fermionic excitations are gapped at all ener- 
gies below A, for collective excitation with energy below 
2A, we can safely integrate out all fermionic degrees of 
freedom and obtain an effective action for the collective 
mode only 







S[9] = dr d 6 x Trln£[0], 



(19) 



where the trace is taken over the fermion momentum, 
frequency and Nambu-Gorkov vector. We have neglected 
here a constant which is irrelevant for the following dis- 
cussions. 

The next task is to expand the action in powers of 9(x). 
We use the standard derivative expansion. With the two 
matrices V\ and V 2 defined in the fermion momentum, 
frequency and Nambu-Gorkov space 0, 
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(20) 



we can expand the effective action to any order of 9. Since 
we consider only the low energy behavior of the collective 
mode, we need only the expansion up to the quadratic 
term, 



S[9] = Tr[GV 2 } + iTr[£Vi£Vi] 



(21) 



After a straightforward calculation, we have 



S[0] = £ 



\0(q)\ 2 (22) 



with the four momentum q = (iw„,q) of the collective 
mode, 77 = n a /m a + iib/nib, and 



V ^ = ^ETr[r 3 «(p + g )], 



(23) 



Now the trace is taken only over the Nambu-Gorkov 
vector. We should note that the functions V and 
V are related to the density and current correlation 
functions [9L fTo| (p(x)p(0)) and (Ji(x)Jj(0)}, respectively. 
Also, we want to emphasize that due to the mass differ- 
ence, the function V is quite different from its usual form 
in the case with equal masses. 

We now discuss the low energy behavior of the collec- 
tive mode. The energy of the collective mode is defined 
through the analytical continuation iu> n — > lo + i0 + ■ In 



the low energy limit, namely ui <§C 2A and «f|p| <C 2A, 
where the Fermi velocity vf is defined as vp = pp/m, 
we can expand the functions T> and V in powers of the 
momentum q. To the leading order, the effective action 
becomes 



S[9] 
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Note that the second term of B comes from 'P(O) which 
vanishes for symmetric systems with a = 1. The energy 
spectrum of the collective mode is obtained by finding 
the zero of the coefficient of the quadratic term in 9, 



(25) 



with the definition for the sound velocity v s = \J B/A. In 
weak coupling limit, the chemical potentials /i a and pib of 
the two species are approximately their Fermi energies, 
and the average chemical potential /x is much larger than 
the gap parameter, the integrations in A and B can then 
be integrated out, 



.4 



mpF 



m 3 



1 



a + 1 



Pf 



(26) 



Together with the total fermion number n — p^/(37r 2 ), 
the sound velocity has the same mass ratio dependence 
as the ratio T c /Aq, 



vf_ ( 2y^ \ 
V3 \l + a) 



(27) 



In the limit a 



— * 1, we recover the well known result 
In the other limit a — » 00, the sound 
velocity tends to be zero. 

At sufficient low temperature T < T c , the specific heat 
c is dominated by the Goldstone mode which has the T 3 
power law 



2tP 
I5vi 
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(28) 



which means that the specific heat for the mixed 6 Li and 
40 K system is about 17 times the one for the pure 6 Li 
system if n is keep fixed. By experimentally measuring 
the specific heat, we can determine the sound velocity 
and check the theoretic result. 

If the fermions carry electric charges, the Goldstone 
mode will disappear due to the long range electromag- 
netic interaction between the fermions [9|. In the lan- 
guage of gauge field theory, the Goldstone mode is eaten 
up by the electromagnetic field via Anderson- tliggs mech- 
anism. This phenomenon is also called Meissner effect. 
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limit a — > oo, the penetration depth approaches infin- 
ity, which means an ideal type-II superconductor. Since 
the penetration depth depends strongly on the mass dif- 
ference, it may change the type of superconductors |l l| . 



FIG. 2: The common a dependence of the ratio of the tran- 
sition temperature to the zero temperature gap, the sound 
velocity, and the inverse penetration depth. 



Suppose the two species carry charges eQ a and eQb, re- 
spectively, the London penetration depth A^ is related to 
the function V\&, 



1 



47TA| 



2m 



(Q 2 a + Ql)+e 2 V(0) 



with the modified matrix 



Qa/m a 
Qb/m b 



(29) 



(30) 



appeared in V due to the charge difference between the 
two species. 

For a neutral Cooper pair with Q a = —Qb = 1, the 
inverse of the London penetration depth squared is zero 
as we expected, 



1 



4^A| 



r dp pi A2 

m 3 ./q 27r 2 m 2 



0. 



(31) 



For the case with Q a — Qb = 1, the mass ratio depen- 
dence is again the same as the ratio T C J Ao and the sound 
velocity v s , 



K 1 



Anne 2 ( 

1 + a 



(32) 



In the limit a — > 1, we recover the familiar penetration 
depth A^ 1 = (47rne 2 /m) 1 ^ 2 . On the other hand, in the 



For systems with fixed reduced mass m (and hence 
fixed vf) but different mass ratio a, from the equations 
(JIJ, 127|) and 132|) . we can define a quantity r which de- 
scribes the common mass ratio dependence of T c /Aq,v s 



and A, 



r c (a)/A (a) = «j(a) = A^(a) = 

T c (l)/A (l) v.(l) A£ x (l) 1 + . 



(33) 



and plot it as a function of a in Fig[2 where 
T c (l), A (l), v s (l) and Al(1) are quantities for the sym- 
metric system with equal masses m a = nib — Tin and 
equal chemical potentials fj, a = Hb = M- 

Based on a general four fermion interaction model, we 
have derived the sound velocity in a superfluid and the 
penetration depth in a superconductor for systems where 
the two species of fermions can possess different masses 
and chemical potentials but keep the same Fermi sur- 
face. We found that the sound velocity and the inverse 
penetration depth have the same mass ratio dependence 
as the ratio of the transition temperature to the zero 
temperature gap. While our result is obtained in weak 
coupling BCS region, we expect that the qualitative ef- 
fect will still work in strong coupling systems, such as the 
mixed atom gas of 6 Li and 40 K. In this Letter, we did not 
consider the Fermi surface mismatch which can induce 
exotic states such as breached pairing and LOFF states. 
In these states, there exist gapless fermionic excitations 
and one can not safely integrate out all fermionic degrees 
of freedom. When we apply the above effective theory 
to such gapless states, we will get imaginary sound ve- 
locity and penetration depth[12j. which have been widely 
discussed in the study of color superconductivity [Hi. \lA\ . 
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